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ABSTRACT 

We use a sample of newly-discovered globular clusters from the Pan-Andromeda Archaeological Survey 
(PAndAS) in combination with previously-catalogued objects to map the spatial distribution of globular clus- 
ters in the M31 halo. At projected radii beyond « 30 kpc, where large coherent stellar streams are readily 
distinguished in the field, there is a striking correlation between these features and the positions of the globular 
clusters. Adopting a simple Monte Carlo approach, we test the significance of this association by computing 
the probability that it could be due to the chance alignment of globular clusters smoothly distributed in the M3 1 
halo. We find the likelihood of this possibility is low, below 1 %, and conclude that the observed spatial coin- 
cidence between globular clusters and multiple tidal debris streams in the outer halo of M3 1 reflects a genuine 
physical association. Our results imply that the majority of the remote globular cluster system of M3 1 has been 
assembled as a consequence of the accretion of cluster-bearing satellite galaxies. This constitutes the most 
direct evidence to date that the outer halo globular cluster populations in some galaxies are largely accreted. 

Subject headings: galaxies: individual (M31) — galaxies: halos — globular clusters: general 



1. INTRODUCTION 

It has long been suspected that a significant fraction of 
the Milky Way's globular clusters formed in smaller 'proto- 
galactic fragments' that were subsequently accreted into the 
Ga lactic potential well. First proposed in the seminal paper 
by ISearle & ZinnI (Il978h . there has since been gradual ac- 
cumulation of indirect evidence in support of this hypothe- 
sis - modern data suggest that the abundances, ages, veloci- 
ties, horizontal-branch morphologies and sizes of many glob- 
ular clusters at Galactocentri c radii > 10 kpc are consistent 
with an externa l origin (e.g.. | Zinnll 19931 : iMackey & Gilmora 
l2004t Marin-Fr anch et afl 120091) . Even so, it has proven 
problematic to unambiguously identify indvidual clusters as 
having been accreted into the Galaxy. The only direct ob- 
servation of this process is the disrupting Sagittarius dwarf 
(llbataet alJll994i) which is depositing at least five globular 
clusters into the Milky Way halo (Da Costa & Armandroff 
Il995t iMartinez-Delgado et al.ll2002i: iBellazzini et al. 2003); 
more controversially, the putative Canis Major dwarf may also 
be responsible for several new arrivals ( iMartin et al.ll2004l) . 
The extent to which the observed properties of sub-groups 
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within the Milky Way globular cluster system reflect the as- 
sembly history of the Galactic halo thus remains a critical un- 
resolved question. 

As the nearest large spiral galaxy, M3 1 is an attractive al- 
ternative target for studying this problem. Its globular clus- 
ter system is several times larger than that of the Milky Way, 
and observational investigation of its halo is less prone to 
the vagaries of projection and extinction that plague Galac- 
tic surveys. It is known that globular clusters projected near 
its central regions exhibit some evidence for sub-clustering 
in position- velocity space that may signal an accretion origin 
(lAshman & Bird 1993; Pe rrett et alJl2003h . although interpre- 
tation is difficult because of the complex nature of the inner 
M31 system. Potentially less confusing are halo regions at 
projected radii Rp > 15 kpc, where dynamical times are also 
longer; however, it is only relatively recently that these remote 
areas have been targeted by deep wide-field surveys. Various 
such studies have shown the M3 1 halo to be littered with co- 
herent tidal debris features indicative of one or more accretion 
events (iFerguson et al ."2002':'Ibata et alJl2007l ). and have also 
facilitated the dis cover y of significant samples of remote M31 
globular clusters jHuxor et al.l l2008) so that it is now possible 
to begin assessing how these objects relate to the stellar halo 
in this galaxy. 

iHuxor et aP (l2010ah have derived the first radial surface- 
density profile for M3 1 globular clusters to extend into the far 
outer halo (Rp sa 100 kpc). Their profile exhibits a distinct 
flattening beyond Rp w 30 kpc, very similar to that observed 
for the metal-poor field halo, and has been interpreted as evi- 
dence that accretion processes have played a role in building 
up both components. 

In this Letter we use new results from the Pan-Andromeda 
Archaeological Survey (PAndAS; McConnachie et al. 200S^ 
to provide the most extensive map to date of globular clus- 
ters in the M31 halo, and explore the implications for the ori- 
gin of this system. PAndAS is an ongoing large program on 
the Canada-France-Hawaii Telescope, utilising the MegaCam 
imager to obtain a deep panoramic view of M31 and M33. 
First-semester imaging and data reduction was completed in 
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mid-2009, revealing in exquisite new detail the abundance of 
low surface brightness sub structure present in the M3 1 halo 
dMcConnachie et al.ll2009i see also Fig. 1). 

2. THE M31 GLOBULAR CLUSTER SAMPLE 

We consider a globular cluster ensemble defined by con- 
firmed obje cts in V3.5 of the Revised Bologna Catalogue 
(RBC; e.g.. iGalleti et al.ir2007lPl plus newly-discovered clus- 
ters from the first semester of PAndAS imaging. 

In this work we are primarily interested in globular clusters 
lying outside Rp = 30 kpc (see below). The RBC list of con- 
firmed globular clusters inclu des 41 objects discovered in our 
pre-P AndAS M31 surveys ( Martin et aT] 120061: iHuxor et alJ 
120081) . of which 3 1 fall beyond 30 kpc. The RBC also contains 
three clusters outside 30 kpc not discovered by us. W e retain 
all RBC entries defined as "extended clusters" (Hux or et alj 
12005, 2008). At present there is little evidence that these are 
anything other than bona fide glo bular clusters with peculiarly 
diffuse structures (IHuxor et al.l 120^1 Oa). at least in term s of 
their constituent st ellar popula tions (Ma ckev et al.1120061) and 
internal dynamics dCollins et a l. 2009). 

The extension provided by first-semester PAndAS imaging 
over the reg i on pre viously surveyed for globular clusters by 
IHuxor et"an ( 120081) consists of nearly complete coverage of 
the inner parts of M3 1 together with a large halo area to the 
west and north-west. We have searched all fields at Rp > 30 
kpc as well as many fields interior to this, using procedures 
similar to those described by Huxor et al. (2008). This has re- 
sulted in a catalogue of 43 previously unknown globular clus- 
ters, the properties of which will be detailed in a forthcoming 
paper (Huxor et al. 2010b, in prep.). For now it is sufficient to 
note that 33 of these objects fall at projected radii beyond 30 
kpc, including 12 between 50-100 kpc and four outside 100 
kpc. Sixteen of our newly-discovered clusters appear to be of 
the extended variety. 

In summary, our sample contains 67 clusters with Rp > 30 
kpc, although only 61 of these lie within the present PAndAS 
footprint. 

All PAndAS imaging is taken during dark sky conditions 
with seeing better than 0.7". Under such circumstances glob- 
ular clusters in the M3 1 halo partially resolve into stars, mean- 
ing that identification is straightforwa rd and unambiguou s 
(see Fig. 1). Following the analvsis of iHuxor et al.l (I2010al) . 
we believe our search procedure does not lead to significant 
bias or incompleteness in the overall cluster selection function 
down to My « -5. Furthermore, away from the very inner- 
most M3 1 fields where crowding is non-negligible we expect 
no significant spatial variation in completeness. 

3. ANALYSIS AND RESULTS 

Figure 1 shows the positions of all globular clusters in our 
sample overlaid on the PAndAS metal-poor ([Fe/H] < -1.4) 
stellar density map. In the outer parts of the M3 1 halo where 
large tidal debris streams are readily distinguished (Rp > 30 
kpc), there is a striking correlation between these features and 
the positions of many globular clusters. Indeed, close inspec- 
tion of Fig. 1 reveals very few remote clusters that do not 
project onto some kind of underlying field overdensity, even 
though these substructures clearly occupy only a relatively 
small fraction of the overall survey footprint. 

' RBC V4.0 was released during preparation of the present work. How- 
ever, tlie updates do not alter our results; in particular, V4.0 does not contain 
any newly confirmed or demoted clusters outside 30 kpc. 



In order to put this result on more quantitative ground we 
undertake several calculations aimed at estimating the proba- 
bility that the apparent association between clusters and debris 
streams could be due to chance alignment. That is, we aim to 
compute the level of significance at which similar substruc- 
tures exist in both the globular cluster and field star distribu- 
tions beyond Rp « 30 kpc. 

3.1. Mock M31 globular cluster systems 

We base our analysis on a simple Monte Carlo method- 
ology utilizing a set of 1.5 x 10^ random realizations of a 
smoothly-distributed M3 1 globular cluster system. These rep- 
resent the null case where globular clusters constitute a well- 
mixed, unstructured halo population. Comparing the proper- 
ties of these mock systems to those of the real M3 1 system 
then allows us to assess whether the observed globular clus- 
ters are indeed spatially correlated with underlying field over- 
densities, or whether the apparent association can be ascribed 
to stochastic effects. 

In each mock system we randomly generated cluster galac- 
tocentric radii between /?p = 30- 130 kpc using a probability 
distribution function defined by the observed globular clus- 
ter radial surface-density profile, and selected position angles 
randomly from a uniform distribution such that each individ- 
ual cluster fell within the PAndAS footprint. Our outer limit 
is the maximum radius with nearly complete coverage over 
the presently-observed area, while our inner limit is the ap- 
proximate radius interi or to which there are many overlap- 
ping halo features (e.g., iFerguson et al.ll2002h and it becomes 
meaningless to associate clusters with field substructure via 
spatial coincidence alone. We set the total number of clusters 
in each mock system equal to the size of the known M3 1 sam- 
ple between /?p = 30- 130 kpc within the PAndAS area (61 
objects). Our derived cl uster surface-density profile closely 
matches that obtained bv lHuxor et al.l (I2010al) . 

3.2. Direct comparison with the stellar density map 

We utilized a FITS version of the PAndAS stellar density 
map to test whether the observed M3 1 globular clusters are 
preferentially projected against regions with higher densities 
of metal -poor red giant stars. The FITS map has an embedded 
World Coordinate System, which allowed us to easily calcu- 
late the mean value in a 7 x 7 pixel box about the position of 
each cluster. The clusters themselves are not generally visible 
on the map since they are usually unresolved by the catalogu- 
ing software. To be sure, we omitted the central (cluster) pixel 
from the average (one pixel k, 350 pc). Our box corresponds 
to a ~ 2.5 X 2.5 kpc region on the sky, representing an ad- 
equate compromise in obtaining enough pixels without con- 
sidering an unreasonably broad area about each cluster It is 
also sufficiently larger than the '--^ 2.5 pixel Gaussian smooth- 
ing kernel used while generating the map. We considered the 
mean value in each box, rather than all pixels individually, be- 
cause this smoothing means that adjacent pixel values are not 
independent. 

We repeated this process for all real and mock globular 
clusters, and formed the results into two cumulative distri- 
butions (Fig. 2). These have quite different shapes - in par- 
ticular, that for the real cluster system clearly contains fewer 
low values than does the distribution for the mock systems. 
In other words, the observed M3 1 globular clusters do prefer- 
entially project onto regions of higher field star density than 
would be expected if they constituted an unstructured halo 
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Figure 1. First-semester PAndAS map of the spatial density of stellar sources possessing luminosities and colours consistent with being metal-poor red-giant 
branch stars ([Fe/H] < —1.4) in the M31 halo I McConnachie et al. 2009). The two dashed circles, representing Rp = 30 and 130 kpc, indicate the vast scale of the 
survey. Our globular cluster sample is overlaid, marked by red points (compact clusters) and blue points (extended clusters). Objects outside the PAndAS area 
ai'e from our previous survey work. Major halo substructres are labelled (see text for details); region (1) indicates the ill-defined major-axis feature and nearby 
overdensities to the east and north, while (2) marks the inner western cluster group. The lower panel shows 1' X 1' PAndAS i'-band thumbnails for ten of our 
globular clusters spanning 30 < Rp < 120 kpc and a wide variety of sizes and luminosities. The lower right-most two are good examples of extended clusters. 
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Figure 2. Cumulative distiibutions of mean pixel values about 61 observed 
M3 1 globular clusters with Rp > 30 Icpc (solid black line) and clusters in the 
1.5 X 10^ mock systems (dashed blue line). The vertical dotted line indicates 
the maximum separation between the two distributions. 

population. A simple Kolmogorov-Smirnov test provides an 
estimate of the significance of this result: the probability that 
the two distributions were drawn from the same parent distri- 
bution is only w 1.3%. 

This straightforward calculation provides a clear quantita- 
tive demonstration that the global spatial coincidence between 
globular clusters and tidal debris streams visible in Fig. 1 is 
almost certainly not due to chance alignment. 

3.3. Stream-by-stream analysis 

Further inspection of Fig. 1 suggests a more complicated 
picture: some streams apparently possess more clusters for 
their size than do others, while several are devoid of clus- 
ters altogether To investigate this in more detail we con- 
sidered each major halo overdensity individually. For now, 
we restrict ourselves to only those stellar substructures pre- 
viously published in the literature. These are highlighted in 
Fig. 1: the gi ant stellar stre am (Ibata et al. 2001); the north- 
east structure (I Zucker et al.ll2004; .Ibata et al..,2005.) : the four 
minor-axis tangent streams A-D jlbata et al.ll2007l) . the latter 
of which extends to a coherent arc to the east of M3 1 ; and the 
southwest cloud and no rthwest stream (McConnachieet al. 
120091) . Ilbata et aP (120071) also describe the "major-axis diffuse 
structure", an ill-defined feature extending to the southwest of 
M31 (region 1 in Fig. 1). Although there is evidently an over- 
density along the major axis here, unlike the other streams 
it is not easily identifiable as a single structure; indeed there 
are additional faint overdensities extending east and north of 
this feature. We defer analysis of this region until it has been 
properly characterized, but note that numerous globular clus- 
ters (w 8 - 10) clearly project onto it. 

We utilized FITS versions of both our metal-poor stel- 
lar density map, and an additional map including stars with 
photometric metallicities -1.4 < [Fe/H] < -0.7, to delin- 
eate the edges of each substructure. These more metal-rich 
stars alter the appearance of several of the o verdensities (e.g., 
Ilbata et al.ll200'7l;lMcConnachie et aO2009h : most notably the 
giant stream, which increases in radial extent and fans west- 
ward; the southwest cloud, which becomes more prominent; 
and stream C, which is considerably broadened - indeed, this 
feature is known to consist of two distinct overlapping com- 
ponents (Chapman et al. 2008). 

For each substructure we measured the mean and standard 
deviation of pixels in numerous nearby regions, and defined 
the edge of the structure by following a contour level w 3.5cr 
above the local background. No globular clusters were over- 
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Figures. Major substructures in the M31 halo. Features associated with 
multiple/zero clusters are outlined in red/blue, while the two cluster overden- 
sities are bounded in green. Clusters are magenta points. The dashed ckcles 
indicate Rp = 30 and 130 kpc. 

plotted during this process. Our results are shown in Fig. 3. 
Although the edges of these substructures are by nature dif- 
ficult to define, they are all sufficiently unambiguous not to 
alter the conclusions we draw below. 

We next counted the fraction of mock systems in which 
A^gc or more globular clusters overlap spatially with a given 
substructure, where A'gc is the number of clusters observed to 
project onto that feature in the real M3 1 halo. We also consid- 
ered the system globally, grouping together all the identified 
substructures. 

Our results are summarized in Table 1 . Taking all substruc- 
tures together, we find it very unlikely that the observed spa- 
tial overlap of 27 clusters with these features can be explained 
by random alignment: the probability sits at just ^ 0.25%. 
This strongly reinforces the result of our previous calculation 
involving the average local stellar densities. 

Individually, the northwest stream and eastern arc are par- 
ticularly well endowed with clusters. The fraction of mock 
systems in which > 6 clusters fall within the northwest stream 
is ^ 2.7%; this falls to below 1% if the additional boundary- 
straddling cluster to the south is included. Similarly, the fre- 
quency with which at least the 1 1 observed clusters fall within 
the eastern arc is ^ 0.5% in the mock systems. 

For the southwest cloud our calculated probability sits at ^ 
2.5%, while for stream C it is less significant at ^ 7.8%. No- 
tably however, stream C is the only case where velocity infor- 
mation exists for both the field substructure and an associated 
cluster (EC4), unambiguously linking the two JCoUins et alj 
l2009l) . 

The giant stream is notable as the only structure where the 
observed number of globular clusters approximately matches 
the number expected in a smoothly-distributed system. Since 
this stream is by far the most luminous stellar substructure 
in the M31 halo, it appears significantly underabundant in 
clusters compared with the streams described above. This is 
perhaps not too surprising: globular clusters located in the 
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Table 1 

Fraction of mock systems matching tlie observed cluster-stream associations. 



Substructure 


Comment 


Wg, 


Number of 
matching systems 


Fraction of 
matching systems 


Global system 

All substructures 


_ 


27+ 


372 


0.00248 



Substructures with multiple clusters 

Northwest stream 

Southwest cloud 
Eastern arc (Stream D) 
Stream C 
Giant stellar stream 



excluding boundary cluster 
including boundary cluster 



6+ 


3999 


7+ 


1101 


3+ 


3742 


11+ 


791 


3+ 


11723 


4+ 


124893 





27383 





138674 





88828 



0.02666 
0.00734 
0.02495 
0.00527 
0.07815 
0.83262 



Substructures with no clusters 

Northeast structure 
Stream A (125 kpc stream) 
Stream B 



0.18255 
0.92449 
0.59219 



Globular cluster overdensities not associated with identified substructure 



Western group 
Northwest group 



any PA; relaxed bounding box 
any PA; relaxed bounding box 



8+ 
5+ 



7033 
8941 



0.04689 
0.05961 



outskirts of the giant stream progenitor may well have been 
stripped away on earlier orbits about M3 1 . Furthermore, the 
globular cluster specific frequencies of dwarf ga laxies span 
a larg e range -- 0- 30 (e.g.. Miller & Lotz 2007; IPeng et all 
120081) . so the progenitor might have possessed comparatively 
few globulars to start with. 

What about the three stellar substructures not associated 
with any globular clusters? Although this occurs commonly 
in the mock systems for all three features, we have already 
assembled ample evidence that clusters are not smoothly dis- 
tributed in the outer parts of M3 1 . Therefore, these three sub- 
structures simply demonstrate that not all field overdensities 
are necessarily associated with clusters. Streams A and B may 
be remains of low-mass satellites, or they may be 'shells' due 
to the impact of a relatively large accreted galaxy (iFardal et alJ 
1^08; Mori et al. 2008). Either way, their intrinsic low lumi- 
nosity likely explains their paucity of clusters. It is also per- 
haps unremarkable that no members of the halo cluster pop- 
ulation are associated with the northeast structure, which ob- 
servations su^gesU^aO^nsi^n^featureJii^ extended 
disk (e.g.. llbata et ani200llRic"hardson et aUlIOOSh . 

3.4. Globular cluster overdensities 

Finally, we highlight two globular cluster overdensities that 
are evident in Figs. 1 and 3 but not obviously coincident 
with underlying substructures. One sits to the west of M31 
at Rj, « 35 kpc, and the other to the northwest at Rp « 105 
kpc. They represent the two groupings with the largest ra- 
tio of local cluster density to the azimuthal average at given 
radius (a factor w 10 enhancement). The grouping with the 
third-highest ratio (« 7.5) overlaps the upper portion of the 
eastern arc. 

We recognize the a posteriori nature of our identification of 
these two overdensities, and attempt to quantify their signifi- 
cance fairly by searching the mock systems for aggregations 
at similar radii but unconstrained azimuth, and with local en- 
hancement ratios > 5. With these relaxed constraints, indica- 
tive probabilities sit at « 5% for both ensembles. 

^ Number of clusters per unit V-band luminosity, normalized at My = —15. 



This is again fully consistent with globular clusters not be- 
ing smoothly distributed in the outer M3 1 halo, although here 
we cannot identify any corresponding field substructures. We 
hypothesize that these two cluster groups may trace under- 
lying substructures that fall below the PAndAS low surface 
brightness limit; confirmation of the nature of these clus- 
ter overdensities will hence require radial velocity measure- 
ments. It is intriguing that the very massive globular clus- 
ter Gl, tho ught to be the stripped core of a former nucleated 
dwarf (e.g.. lBekki & Chibal2004^ . is a prominent member of 
the western group. 

4. DISCUSSION 

Put together, our results provide strong evidence that glob- 
ular clusters in the outer halo of M3 1 are anisotropically spa- 
tially distributed, and preferentially associated with underly- 
ing tidal debris features. Of the 61 clusters in the PAndAS 
footprint with Rp > 30 kpc, a t lea st 27 lie on the major sub- 
structures outlined in Section [331 if the complex major-axis 
region is also included, this number rises to « 37. A fur- 
ther 13 objects are memb ers of the two cluster overdensities 
described in Section [341 leaving just 11 of the sample unac- 
counted for. Figure 1 reveals only a handful of these lie away 
from field overdensities altogether 

These numbers imply that the majority (> 80%) of the 
outer globular cluster system of M3 1 has been built up via 
the accretion of satelUte host galaxies. This fraction matches 
closely that inferred for the outer Galactic system (e.g., 
iMackey & Giknorel20 04' Forb es & Bridg es 2010). Ou r work 
provides a striking direct illustration of the Sea rle & ZinnI 
(119781) paradigm; along with the disrupting Sagittarius dwarf, 
it represents the most clear-cut observation to date of the as- 
sembly of a globular cluster system in action, and does so on 
a grand scale - across an entire galactic system and a huge 
population of clusters. 

We note that all but one of the extended clusters in our sam- 
ple with Rp > 30 kpc are projected onto stellar substructure or 
are members of a cluster overdensity. This observation is con- 
sistent with the idea that these puzzling objects may predom- 
inantly originate in lower-mass galaxies (e.g.. .Da Costa et aJl 
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l2009h . 

Intriguingly, streams without clusters are seemingly in the 
minority in the outer M3 1 halo. With the caveat that there may 
be very low-luminosity features falling below the PAndAS 
surface-brightness limit, this could imply that most substruc- 
ture in the outer M31 halo is due to the accretion of just a few 
larger cluster-bearing satellites. This would be in qualitative 
agreement with the results of cosmologically-motivated sim- 
ulations, which suggest that the haloes of large spiral galax- 
ies are predominantly assembled via the accretion of a hand- 
ful ( < 5) of significant pr ogenitors (e.g., Bullock & Johnston 
IM)llCoopCTerani2M0h . 

Owing to their recent discovery, little is yet known about 
most of the remote M31 clusters described in this pa- 
per However, a few have published colour-magnitude di- 
agrams and/or spectra, which reveal them as almost exclu- 
sively metal- poor with —2.3 < [ Fe/H] < -1.6 (M ackey et al.l 
12006. 20071 120101: lAlves-Britoet al. 2009). Others exhibit 
a small dispe rsion in integrated colour, further supporting 
this assertion (^ Huxor et al.ll2010ah . This might imply that a 
significant fraction of metal-poor globular clusters in large 
galaxies are accreted objects, as suggested by a variety of 
models (e.g.,|Cpte et al. 1998, 2000rPrieto & Gnedin 2008; 
iMuratov & Gnedi32010l) . Our present results raise the excit- 
ing imminent prospect of characterizing, for the first time, the 
individual properties of multiple accreted families of globular 
clusters in a galactic halo. 
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